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Abstract 

\& , We consider the stationary incompressible Navier Stokes equation in the exterior of a disk Bel 2 

with non-zero Dirichlet boundary conditions on the disk and zero boundary conditions at infinity. 
We prove the existence of solutions for an open set of boundary conditions without symmetry. 
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1 Introduction 

In this paper we consider the incompressible Navier Stokes equations in an exterior domain: 



Au-Vp^u-Vu, inR2 ^ (1) 
div u = , 

with B a smooth bounded domain, with non-zero Dirichlet boundary conditions on dB, and zero boundary 
conditions at infinity: 

u !aB = u* , lim u(x) = . (2) 

|x|— ►oo 

Of particular interest is the case of boundary data u* with zero flux: 

u* • n da = . (3) 

dB 
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We note that, since the size of B is arbitrary, we have set without restriction of generality all the physical 
constants in ([TJ) equal to one. 

The above system is a special case of the exterior Navier Stokes problem: 

| -(u.V)u-A9 lU + Au^Vp = 0, inR n N5j (4) 

with n — 2 or 3, with B a smooth bounded domain in W l , with boundary conditions @, and with 
A G {0, 1} distinguishing between the case of a flow "around" B (A = 0) and a flow "past" B (A = 1), 
respectively. The system ([IJ-© corresponds to n = 2 and A = 0. The case A = is in many respects 
more complicated than the case A = 1, and, whereas the picture is rather complete for n = 3, the case 
n = 2, A = 0, presents particular difficulties. The difficulty with the classical method for solving the 
Navier Stokes equations consists in the fact that the linearization around u = is given by the Stokes 
system, which, for n = 2, does not admit a solution satisfying ([2]), unless the domain B and the boundary 
data u* satisfy certain symmetry conditions. This fact is known as the Stokes paradox. For completeness 
we note that if one relaxes the no flux condition @, there exists a two parameter family of solutions to 
([I])-©, the so called Hamcl solutions, see [7]. These examples emphasize that the decay of solutions can 
be arbitrary slow and that uniqueness might be lost for some boundary data. However, these solutions 
have flux larger than one, and are far from the regime which we will consider here. 

In what follows we construct a new class of solutions to (HJ-©, by linearizing not around u = 0, but 
around u = fix 1 - / |x| 2 , with > a/48. This improves the decay of the solutions to the vorticity equation, 
yielding vorticities decaying at infinity generically faster than |x|~ , instead of like |x|~ as would be the 
case for the Stokes equation, thus avoiding the Stokes paradox when reconstructing u via the Bio-Savart 
law. 

To put our problem into a wider context, we briefly recall the concept of weak solutions for (|1]), $2$ 
(also known as generalized solutions or D-solutions), and the method of J. Leray [TB] for proving the 
existence of such weak solutions. 

Definition 1 Given u* £ H 1 / 2 (dB) satisfying (0], a function u which satisfies the following conditions 
is called a weak solution to (0): 

1. u G D 1 ' 2 (R n \B), where D 1 ' 2 (R n \B) is the subset of Lj oc (R 2 \B) containing functions with gradient 
in L 2 (R n \B), 

2. u is divergence-free and u = u* on dB, 

3. for all divergence- free vector fields w € C^°(R n \ B), there holds: 

/ _Vu:Vw+/ _((u ■ V)u + Adiu) • w = . 

JR"\B JR"\B 

The method of J. Leray to prove the existence of solutions according to this definition, and a posteriori 
to ((4|), ([2|), in the sense of distributions, consists in the following steps: 

• First, one introduces a sequence of approximate problems by restricting Q to bounded subsets 

C M" containing B, with zero Dirichlet boundary conditions on dil \ dB. 

• Second, one proves the existence of (weak) solutions to all these approximate problems. 

• Third, one shows that for any sequence of bounded subsets exhausting W 1 \ B, there exists a 
subsequence, such that the corresponding approximate solutions converge to a weak solution of (j4]) , 
©• 

• Finally, given a weak solution u, a pressure p an be constructed via De Rham's theory, such that 
the equations ^ are satisfied in 2?'(R n \ B). 
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See also [T2J H31 HOI H31 HI] j where this method has been adapted to a similar system with more general 
boundary conditions. Note that if B has a smooth boundary, the ellipticity of the Stokes operator (see 
[5j Section IX. 1]) and the smoothness of u* imply that weak solutions are smooth. Therefore, for smooth 
data, the only possible shortcoming of weak solutions is that they may not satisfy the boundary condition 
at infinity in a point-wise sense. Much work has been devoted to clarify the situation in various cases 
(see [7] for more details): 

For n = 3, the condition u £ D 1,2 (]R 3 \ B) implies that weak solutions tend to zero at infinity. The 
exact decay can be obtained by various methods yielding the following results: 

• for A = 1, there exists a solution that decays like the fundamental solution of the Oseen equation 
(the linear system obtained from ((4]) by deleting the nonlinear convective terms) [21 [3j [4] . This 
result can be obtain by a detailed analysis of the Oseen equation with a source term in the usual 
Sobolcv spaces [U 0] , and also in weighted Sobolev spaces [3] . 

• for A = and sufficiently small boundary data, there exists a unique weak solution, and this 
solution decays like a Landau solution |15) , a special solution of the nonlinear system which decays 
like l/|x|. This result can been obtained by constructing first a strong solution to ([2]), (jU, which 
is asymptotic to the Landau solution, by perturbativc techniques. Using the known decay of this 
particular solution as an input [7J Section IX. 9], one then proves a weak-strong uniqueness result 
for small data. 

For n = 2, the situation is more delicate since the condition u € D 1,2 (R 2 \ B) does not guarantee that 
the boundary condition at infinity is satisfied: 

• For A = 1, the relevant linear system is again the Oseen equation, but the results concerning the 
decay are limited to small data, since, as for the case n = 3, A = 0, perturbative techniques are 
used to prove the existence of a strong solution decaying at infinity like the fundamental solution 
of the Oseen equation. This solution is then again used as an input to a weak strong uniqueness 
argument in order to show the decay of weak solutions. These results can be found in [6]. 

• The case A = remains largely open. As we already pointed out, the problem is that the solution 
to the Stokes equation with boundary data u* 7^ diverges at infinity, unless one makes additional 
assumptions on the domain B and the data u*. Partial results for the Navicr Stokes system with 
symmetric data can be found in [8j [18j [19j [XT] ■ 



From now on we limit the discussion to the case where B is a disk of radius one. We choose x = (x, y) 
Cartesian coordinates with the origin at the center of B, (r,8) g Q := (0, 00) x (— 7r,7r) the associated 
polar coordinates, and (e r ,eg) the corresponding local orthonormal basis. For the function u we have in 
polar coordinates: 

u(r,0) = u r {r,6)e r + ue(r,6)e e , V(r,0)efi. (5) 
The following theorem is our main result: 



Theorem 2 Let /io > Merit = v48 and u* <G C°°(dB) satisfying ([3]) be sufficiently close to u* o :=fj,Qeg. 
Then, the equations {!]), H|), with boundary condition u* , have at least one solution {u,p) € C°°(M 2 \ 
B) 2 x C°°(R 2 \ B). Moreover, there exist /1 close to jiQ such that: 



lim 



n{r, 6) -—\L°° =0. (6) 



Remark 3 If the pair (u(x,y),v(x,y)) is a solution for the boundary condition (u*(x,y),v*(x,y)), then 
the pair (u(x, — y), — v(x, — y)) is a solution for the boundary condition (u*(x,~y),~v*(x,—y)). Thus, 
our result extends to fig < ~fJ-crit- 

Remark 4 Ifu(r,9) is a solution for the boundary condition u*on the complement of the unit disk, then 
for all A > 0, Au(Ar, 9) is a solution for the boundary condition Au* on the complement of the disk of 
radius A -1 . 
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Remark 5 The restriction to the case where B is a disk is for the sake of simplicity only. This permits 
to rewrite the system in polar coordinates, yielding explicit expressions for the solutions. We expect that 
with more work the results can be generalized to arbitrary smooth B. 

To prove Theorem [5] we proceed as follows: We fix pL > fi cr u and consider the pair (u^,^): 

u M (r,0) = — , P„(r,0) = ~^, V(r,0)efi, (7) 
r zr 

which is an exact solution to (JTJ) , ([2]). Next we set, (u,p) = (u M + v,p = p^+ q) and prove, that for all 
sufficiently small boundary conditions v* satisfying 

/ v* • n dcr = , (8) 

JOB 

there existence of a solution (v, q) £ C°°(R 2 \ B) 2 x C°°(tt 2 \ B) such that 

v lss — v * + M*; f° r some 

/U* > pL cr it depending on \i and v*. In a final step, wc show that this function can be inverted, giving fi 
as a function of /it* and v* , thus yielding Theorem [5J 

The feasibility of our approach relies on the fact that the system obtained by linearizing (fTJ), @ around 
the explicit solution (u^p^) can be analyzed explicitly. As mentioned above, when compared with the 
case [i = 0, i.e., the Stokes equation, the vorticity decays for p. > fJ, C rit faster than 1/r 2 , instead of like 
1/r, such that u can be shown to decay faster than 1/r at infinity, making the nonlinearity subcritical. 
Introducing suitable function spaces, we are then able to solve the full non-linear system by a classical 
fixed-point argument. 



2 Dynamical system formulation 

Let (u,p) £ C°°(1R 2 \ B) be a solution to |T]), satisfying ([3]). We first make the construction of 
the stream- function if) associated with u precise. Since u £ C°°(R 2 \ B), we have in particular that 
u* £ C°°(dB). Let Mint £ C°°(B) satisfy Ui nt = u* on dB. Such a function exists since u* satisfies ©. 
For instance, Mint can be the solution to the Stokes equations on B, with boundary condition u* on dB. 
Then, setting: 

u in R 2 \ B~, 
Ui n t in B , 



u(x,y) 



we obtain a continuous divergence-free vector-field on the whole of R 2 . Furthermore, this function is 
smooth on both sides of dB so that there exists V £ C^M 2 ) n C°°(B) n C°°(M. 2 \ B) satisfying u = 

Instead of (jT|), ([2]), wc consider now the equation for the stream function ip and the vorticity cj = V x u, 



in M. \B . 

For the function u wc have in polar coordinates (|S|), and the vorticity becomes: 

ui = —d r (rug) dgu r , V (r, 0) £ fi . 

r r 

For the boundary data we have: 

u*(0) =u* r {6)e r + u* e {6)e 6 , V 6 £ (-tt.tt). 
In polar coordinates we get the following equations for the stream function ip and the vorticity to: 

d rr lp + -d r l}j + \dgelp = -0J , 

1 1 u e V(r,0)6fl, (9) 

d rr oj H — 9 r o; H — ^deeuj = u r d r oj H 9ew , 
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and 

Ur r ' V(r,ff)efi, (10) 

ue = -d r if> , 

together with the boundary conditions: 

f u r (l,9) = u*(0) , Hindoo u r (r, 9) = 0, 

< VS6 (-7T,7T). (11) 

[ tte(l,0) = ^(0), Um r ^.oo ue(r, 0) = 0, 

For the exact solution (u M ,p M ) given by ||2J) we have in polar coordinates for the corresponding stream- 
function-vorticity pair (^ M ,w M ), for all fi € K: 

[ <*W,0) = , 

In order to prove Theorem [5] we construct, as explained above, a solution which is a perturbation of the 
explicit solutions (u^,^). We therefore set ip = ip^ + 7 and uj = lu^ + w. Substituting this Ansatz into 
©1 POP - we obtain the following equivalent system for the unknowns (7,u>): 

drrl + 7<9 r 7 + = -w , 

\ fi d el d rl VM)efi, (12) 

d rr w + ^d r w + -irdegw ^dgw = d r w —dew , 

with the boundary conditions: 

3 e 7(M) = v* r (0), 

dr7(l,0) = -v* e (e), V0G(-7T,7r), (13) 

Um^oo (\ 7 (r,e)\ + \dMr,e)\) = 0, 

for certain (v*(0),Vg(0)) to be defined later on, satisfying: 

v* r da = , (14) 

OB 

and which are small in a sense to be made precise. 

Following the method developed in [TJ], we solve (fT2"j) . (fT5|) . for data (w*,^), by interpreting the 
radial coordinate r as a time and by expanding in a Fourier series: 



7 (r, 0) = £ 7 „ (r)e ine , u,(r, 0) = £ W „(r)e 4 ' 



Notation. To unburden the notation we write for the Fourier series of 7 and w: 

7 = ( 7 ™)™ez, u> = (w„)„ e z, (15) 
and analogously for all other functions. 

From (fT2|) . (p~3|) we obtain, for n £ Z, the following system of ordinary differential equations: 

1 n 2 

<9rr7n + -^r7« o7n = ~W n , 

r ' 2 on (l,oo), (16) 

_ 1 _ iiin + n 

O rr Wn + -O r W n IU„ = /<„ , 



with the source term F n given by: 



F„ = — ^2 {kwkd r ji - ljld r Wk) , (17) 

k-\-l—n 
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and with the boundary conditions: 



iwy n (l) = u*„ , 

-^7n(l) = v* e>n , VneZ\{0}. (18) 

lirn^oo (|7„(r)| + \d r j n (r)\) = , 

Note that v* = by assumption p4[) and that the value of 7o(l) is irrelevant, i.e., the stream function 
is only unique up to an additive constant. As we show later in this section, the value v@ cannot be 
chosen freely if one wants the solution 70 to satisfy the boundary condition at infinity. 

For convenience, we first solve (|16[) with boundary conditions: 

7n(l) = In , 

< w n (l) = w*, VneZ\{0}, (19) 

linv^oo (|7 n (r)| + \w n {r)\) = 0, 

instead of (|19p . Once the solution is constructed we then re-express the solution in terms of the original 
boundary conditions. 

Assuming that the functions F n are continuous and decay sufficiently rapidly at infinity, there exits 
exactly one solution to §TE\i satisfying ([TO)) . Since the Green's function of equations p^|) are r 1— > r ± l"l 
and r 1— > r ±< >™, respectively, where ~ y n 2 + i/j/n, with 7Ze(y/z) > for z G C \ (—00, 0], the solutions 
are given by the following explicit expressions: 



with: 




For n ~ 0, there still exist solutions to (|T5)) decaying at infinity, but these solutions exist only for exactly 
one boundary condition. The reason is that for n — the Green's functions for the equations in (fT()]) arc 
r h- > 1 and r t— > lnr, which do not decay at infinity. The solutions decaying at infinity are: 

7o (r) 

< 

w (r) 

We recall that the value of 70(1) is irrelevant. The value of Wq(1) fixes the value of d r 70(1) = 

function of /i. Once the solution is constructed, we will show that /Lt i— > fx — d r "/o(l) =: fi* can be inverted, 

which then shows the existence of a solution for an open set of boundary conditions. 

3 Functional framework and main result 

We now introduce the function spaces which we use to solve the system (TT8"]) -(f2^ j) . We use the notation 
introduced in (fT5l) : 



oc -1 /-oo 



t u>o (t) dt ds , 



00 -1 poo 



tF (t)dt ds. 



(22) 
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Definition 6 Given n > 0, a > and m G N, such that m < k, we set: 

B K := {ip* e C z such that sup„ eZ (l + |n|) K |y>*| < oo} , 
:= {ip* £ B K such that tp^ = 0} , 

and 

B a>K := {ip G (C([l,oo);C)) z , such that sup„ eZ sup re[loo) r Q (l + |n|) K |</?„(r)| < oo} , 
:= G (C m ([l,oo);C)) z , svcA iftai (d l r <p n ) neZ e /or all < I < m} . 

These function spaces are reminiscent of weighted Sobolev spaces, and permit to obtain sharp estimates 
on the decay of solutions to (|20 p -(|22 p . The spaces with one lower index (mainly B° K ) are used for the 
boundary data, whereas the spaces with two lower indices (mainly U™ K ) will be used for solving (f2TJ]) -([2"2 j) . 

The spaces introduced in Definition [6] satisfy the following straightforward properties. Given a > 0, 
k > 0, and m € N, such that m < k, we have: 

1. The spaces B K , B a ^ Kl U™ K are Banach spaces when equipped with their respective norms: 

\\w ■ B K \\ = sup(l + M) K K|, \\0 : B a , K \\ = sup sup r a (l + \n\) K \ Vn (r)\ , 

neM nGZr€[l,oo) 
m 

2. Given a > a' and k > k' we have the embedding S QjK C B a ' , K > together with the bound: 

\\$;B a , jKl \\<\\y;B a , K \\, V£e£ Q ,«. (23) 

3. The space B® K is a closed subspace of and thus also a Banach space. 

We now formulate the problem of finding a solution to (TTg)) - (f2"2")) in such a way that we can apply the 
inverse map theorem on our function spaces: 

Lemma 7 Given /i > Hcrit; let a > be sufficiently small and k > 0. Then, the map S^: ^4+2 Q ,k X 
B® + 4 x B® +2 — > K+4 x W^ +2 K+ 2, which associates to the triple (F, 7*, u)*) £/ie pair (7,10) virtue of 
equations i20\) , (2ty) , together with H9\) , \21} ), is linear and continuous. 

The notion of a small enough will be made precise in the last section. 
Lemma 8 Let a > and k > 0. Then, the map AfC: (U^ K+4 x U^ +2 K +2) 2 ~~ * B±+2u.k+i, defined by 

M£[(%,w a ),(%,w b )}= L^-~J2 {kw a k (r)d rl t(r)~l^(r)d r w b k (r))) , 

\ k+l=n J neZ 

is bilinear and continuous. 

The proofs of these lemmas are postponed to Section [5] We also introduce the trace operator Ti: 

Ti: C X C &Ki-i x &Ki-i j 

(7, w) I — S> ((m7„(l))„ 6Z , ((S n . Q - l)(9 r 7n( 1 ))«6z) , 

where 5 n>m is the Kronecker symbol. This map is linear and continuous for arbitrary (cti,Ki) £ (0, oo) 2 , 
(i = 1,2)'. 
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To compute solutions to (|T51) - (j2"2")) . we introduce a map which allows to solve the differential 
equations and constrain the trace on r = 1 in one step. Namely, given [i > [i C rit, ot sufficiently small and 
k > 0, we set: 

<TV : « K+4 x Ul +W ) x (B° K+4 x B° K+2 ) — ► (Z< K+4 x Z^ +w ) X (B° K+3 x S° +3 ) 

i=(7,t&) \ /5 M (AA£(i,f),x*)-£A ( 24 ) 

x* = (r.«>*)J ^ \Ti[S lt (ArC(&,£),£*)]J 

By definition, if (x = (j,w),x* = (7*, a/)*)) is a solution to $^(3:, x*) = (Q,(v*,Vg)), then (7,10) satisfies 
(fTB|) - ([T5)) . This motivates the following notion of K-solutions: 

Definition 9 Given an exponent k > 0, an angular velocity fi > fJ, C riti o-nd a boundary condition v* := 
(v*,ig) G &K+3 x we ca M K-solution for the boundary condition v* and the asymptotic angular 

velocity \i a pair x = (7, w), such that, for sufficiently small a > and some x* = (7*, w*) G £>„+4 x ^k+2 : 

• X G K a+4 , )K X Wq +2jK+ 2: 

. * M (£,i*) = (0,v*). 

The remaining sections are devoted to the proof of the following result: 

Theorem 10 Given k > and /xo > £t C ri« there exists e K ^ > and an open interval I Ktfi0 3 fJ-a such 
that, given v* G B(B° +3 x $° , 3 ; e KjA10 ) and \l* G I K .fi , £/iere exists /1 > Hcrit an d a K-solution (7, w) for 
the boundary condition v* and the asymptotic angular velocity /1, satisfying the condition fi—d r "fo(l) = (i* . 

As mentioned above, the notion of a small enough will be made precise in the last section. Before 
entering into the details of the proof of Theorem 1101 we explain why it implies Theorem [5] 

Proof of Theorem^ Let (1q > fi cr it and k > 1. Applying Theorem 1 1 01 yields a ball of initial conditions 
with positive radius e KjAlQ > and an open neighborhood I K ,n of (Iq. 

For u* G C°°(dB), we define : 

1 f 2n 
M = it / "e0)ds, 
2ti" Jo 

and the sequences w* and u* by uq = vq = and: 

u*J6)e- ind d9 , v* n = — u e {8)e- m9 d8, VnGZ\{0}. 

2tt Jo 

The regularity of u* yields that (u*,v*) G £>° +3 x £?° +3 and 

u* r {6) = <e in6 , u* e (8) = /i* + <e ine ■ 

We now assume that: 

|| («*,»•) ; £° +3 x S° +3 || < £ft)M , /x* G J« )W) , (25) 

which makes the meaning of u* sufficiently close to fioeg in the statement of Theorem [5] precise. 

Consequently, the assumptions of Theorem [TU] are satisfied, which yields that there exists fi > fx cr it 
and a K-solution (7, w) for the boundary data (u*, v*) satisfying fi — 9r7o(l) = A**- Let 

u;(r,0) = 5>„(r)e <nfl , 7 (r, 0) = ]T 7 n(r)e in9 , V(r,0)efi. 

Because k > 1, classical results from the theory of Fourier series yield that: 

• we C 2 (R 2 \B) and 7 G C 4 (R 2 \S) so that v = G C 3 (R 2 \B), 
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v- Vw G C^R 2 ^) with: 



v • Vw(r, 0) = 



Let 



-- ^ (lwi(r)d r jk{r) - k-f k (r)d r wi(r)) 



k+l=n 



e me , V (r, 0) G Q . 



Aw--e 9 -Vw-v-Vra=:^£ C([l, oo); C([-tt, tt])) . 



Because (7,w) satisfies (fTB")) , all the Fourier coefficients of ip vanish identically on [l,oo). Hence, (7, if) 
is a solution of 

{A7 = — w, 
a H Y7 T7 inK 2 \B. (26) 
Aw eg ■ Vw — v • vw , v v ' 

r 

Straightforward manipulations of the Fourier series of 7 yield that: 

h ; ul, K+i \\ 



|v(r,0)| < 



ra+l 



, V(r,0)efi. 



(27) 



Since u* has zero flux, i.e., since has zero average, we have that v(l,(9) = v*(6)e r + Vg(6)e$, for all 
9 G (— 7r, 7r), where: 



<(0) = £m7„(l) e irie = E< emS = U *M > 

nGZ nSZ 

0) = - £ 9r7n(l)e me = - a r70 (l) - M* 



Therefore, if we set ip := Vv + 7, w := w, u := V ^ = u m + v j then the pair (ip,co) is a solution to 



Aip = -oj. 



in R 2 \ B , 



(28) 



Acj = u • Vw , 

and the following boundary conditions are satisfied (recall that n — <9 r 7o(l) = /U* by construction of 7) 

u = u* , on dB , lim \u(r, 0)| = . 



The inequality (|27p implies that the boundary condition at infinity is satisfied in the following more 
precise sense: 



lim 

r— ¥00 



r (u(r, 6) 



;i°°((-7r,7r)) = lim r\\v(r, 0); L°°((-tt, tt)) || = . 



To complete the proof, we need to show how to obtain the Navier Stokes equations ([T]) from the 
relations between u, ip and uj, together with (|2"5]l . First, multiplying (|2"5]l by <p G C^°(R 2 \ £?) yields: 



2 \B 



2 \B 



[u • Vw] (p . 



We have the basic identities: 

w = V x u , u • Vw = V x [u • Vu] , 
from which we obtain, after integration by parts, that for any given ip G C^°(R 2 \ B) 

- [ _Vu:VVV= / _[u-Vu].V^. 



(29) 



(30) 



This identity yields the pressure p via De Rham's theory, modulo the difficulty, that not all the divergence- 
free velocity-fields w G C^°(R 2 \ B) of compact support can be written in the form V with <p G 
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C^°(R 2 \ B). More precisely, if a smooth velocity-field v satisfies V x v = in R 2 \ B, then v is the 
gradient of a function up to a contribution of the form Cx. 1 - /\x.\ 2 . We now show that this contribution 
vanishes in our case. 

Let $ € C°°(R) be such that supp($ ) CC (1,2) and $ (0) = 0, $ (2) = 1, and let w = V i $ - Jn 
polar coordinates we have w (r, 9) = — & Q (r)eg, for all (r, 0) £ il. Given a divergence-free w £ C^°(R 2 \B) 
we define w and tp by: 



w = w — 7\/ w wq with 



wg(r, 0) dr 



tp(r,6) = / u;9(«,ff)da-M w *o(r). 



By definition of M w , we have that £ C^°(]R 2 \ B) and Y/- 1 ^ = w, so that (f3"0"j) implies: 

- / _ Vu: Vw = / _ [u ■ Vu] • w . 

Jr 2 \b Jr 2 \b 

Replacing w by its definition yields that, for any divergence-free w € C£°(]R 2 \B), we have: 
-/ _Vu:Vw=/ _ [u • Vu] • w-M w / _ ([u • Vu] • w + Vu: Vw ) . 

JR 2 \B JR 2 \B JR 2 \B 

Let Iq be the last integral in the previous equality. We then have: 

Iq = lim / ([u • Vu] • w + Vu: Vw ) . 

N ^°° Jb(R 2 ,N)\B 

Integrating by parts, we obtain, for all N > 1: 

([u • Vu] • w + Vu: Vw ) = / ([$ u • Vu] • n ± +$' d r u ■ n- 1 ) da 

JdB(R 2 ,N) 

- / ($ u ■ Vw + Vw ■ V$ ) 

JB(R 2 ,N) 



B(R 2 ,N)\B 



[ ($ [(u • Vu) ■ n- 1 - <9 n w] +$ <9,.u ■ n x ) da , 

JdB(R 2 .N) 



'dB(R 2 ,N) 

where, in order to get the last identity, we have again used that u • Vw = Aw, in R 2 \ B. Since u decays 
like 1/r, Vu decays like 1/r 2 , and Vw like 1/r 3 . This yields that Jo = in the limit N —> oo. Finally, 
we have: 



/ _Vu: Vw= / _[u-Vu]-w, 

Jr 2 \b Jr 2 \b 



2 \B JR 2 \1 

for any divergence-free vector-field w £ C C °°(IR 2 \B), and Dc Rham's theory (see [2TJ Remark 1.5]) implies 
the existence of a pressure p such that (JTJ) is satisfied. ■ 



4 Proof of Theorem [10 



In this section k > and /io > ^ c nt are fixed. First, we set /i_ = (/.to + Hcrit) /2 and //+ = (2/io+ /i crj f)/2 
so that / := [/i„,/i + ] satisfies 

M0 G [/*->/■*+] C (Merit, OO) • 

We also set: 

lV2 \ 



1 . 

a := — mm 
4 



1 



^l + lM-l' + l 



-2,1 



We emphasize that, because hq > fi C rit, we have a £ (0, 1/4]. Let: 



4>, 



(^a,K+4 X ^a+2,K+2) X X ^ (^a,K+4 X ^a+2,K+2) X (^k+3 X ^"+3) 



f S^ l {N'C{x,x),x*) - x s 



(31) 



(32) 
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We will show in the next section that for [i G I the map <£> M is well defined. We split the proof of 
Theorem [TU] into two steps. First, we show that we can construct a K-solution for any sufficiently small 
boundary condition (u*,v*) G <6«+3 x @k+3 an d an interval of asymptotic angular velocities /i. In a 
second step, we analyze the dependence of this solution as a function of \x. 

We have the following abstract result: 

Proposition 11 Let 1 be a compact interval and X , Y two Banach spaces. Assume that $ : fx G 1 1— > 3"^ 

satisfies: 

. ^^Cil-C^X-Y)), 

• = for all fi £ I, 

• _D$ M (0) is one-to-one and onto, and has a continuous inverse for all /i G /. 

Then, there exist positive constants r\ x and r\ y , such that $^ is a C 1 -diffeomorphism from Bx(0, r) x ) onto 
*JV (^-^(0, f]x)) -> By(0, rjy). Furthermore, the family of inverse maps $ _1 : fj, G 1 1— > ( < i ) A ,) 1 satisfies 

Proof. The proof is standard, but for the sake of completeness we recall the main ingredients. Given 
the hypothesis of Proposition [11] the map $ M satisfies the assumptions of the inverse function theorem 
for arbitrary fj, G I, so that there exits n Xtfl > and rj y ^ > such that $ M is a C 1 -diffeomorphism 
from Bx(0,ri Xyfl ) onto <S>^ (Bx(0, r] Xtli )) D By (0, r] M ). Since <I> is continuous, it is clear that these 
constants can be chosen independently of /x, locally in fi. By a compactness argument, we can therefore 
find constants r\ x > and ij y > such that $ M : Bx(0,ri x ) — > ^ (Bx(0, n x )) D £?y(0,n y ) is a C 1 
diffeomorphism for arbitrary /i G /. We now show that S" 1 G C(I x ^(0,%)). The proof that 
G C(/; C 1 (By (0, 77j,); -Bx(0, ?7a;))) is then obtained by differentiating (with respect to x) the identity 
o ^^(a;) = x which holds true on Bx(0,n x ). 

Given (fi, ft) G I 2 and {y,y) G [-By(0, r? y )] 2 , we denote: 

By construction we have: 

y - y = ^0) - <&fi{x), 

= D$„(0)[ai -x] + - * A (5) + ofla: - £|) 

= D$ M (0)[z - 5] + o(| M - + o(||as - i; X||) . 

Consequently, reducing the size of n x and rj y if necessary, we get that: 

\\x-x;X\\ < 2\\[D^,(0)]- 1 :C C (Y;X)\\ [\\y - y; Y\\ + o(|/i - £|)] , 

where £ C (F; X) denotes the set of continuous linear map Y — > X. This completes the proof. ■ 

We now show that we can apply Proposition QT] to the map as defined in To this end, we 
remark that $ M depends on /i only through \i i— >■ 5^, and that for all /x G / the differential £>$ M (0) is: 

D%(0)[x,x*] = (1^(0, **),-*), V(x,F) G « K+4 x^ +2iK+2 ) x (6° +4 x S° +2 ) . 

Let Tfj,(x*) := rnS M (0,x*). Since $ is a combination of linear and bilinear maps, it suffices to apply 
Lemma [5] and the following lemma in order to check that the assumptions of Proposition 1 1 1 1 are satisfied: 

Lemma 12 Let a be given by (|3ip . then the restriction of the map S: /i <— > to I := satisfies: 

i) Se C(L;£ c (B K , 2a+4 x (B° K+4 x B° K+2 ) ; < K+4 x Z£ +2)B+2 )), 

7^ is a one to one and onto map B® +i x £^ +2 — ^ x ™^ continuous inverse. 
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We postpone the proof of this technical lemma to the next section. We now apply Proposition [TT] to 
$, but restrict the image to the component x. This yields the following result: 

Theorem 13 There exists a map \&: fx i— > "I^ satisfying 

* 6 C(I;C\B B o +3XB o +3 (0,e K ) ; < K+4 x U 2 a+W )) , 

such that, for all v* G £? B o 3 xB° 3 (0> e «) fflnc ^ M € /, ^(v*) is a K-solution for the boundary condition 
v* wzi/i respect to the asymptotic angular velocity /i. 

In a final step we show how to prescribe the zero mode of the solution ^ (v* ) by using the dependence 
of the solution on \i. Let 77 = (p.Q— At cr i*)/4, and consider a boundary condition 

Let (VVjWju) := v E' M (v*), for all /j, E I. Using that ^(O) = 0, and restricting the size of e K if necessary, 
we can assume that 

|d r ^,o(i)|<||Vv;<J|<??, Vpei. 

Consequently, the map fi h-> /i — d r V7j,o(0), which is continuous from / to R, because fi H> \& M (0, v*) is 
continuous, satisfies: 

/i_ - ar^_,o(l) < M- + »7 < Mo , M+ - 9r^+.o(°) > M+ - ^ > A*0 ■ 

Hence the image of this map contains an open interval I K ,it containing ji®. This completes the proof of 
Theorem [TDJ 



5 Proof of main lemmas 

This section contains the proof of the technical lemmas which have been used without proof in the 
previous sections. First we prove Lemma |S1 which is standard. We then give proofs of Lemmas 171 and [T2l 
which are more delicate. 

Proof of Lemma\E Let F = Af£{(j a 1 w a ),{j b ,w b )) 1 for (f ,w l ) € U 2 a K+A x U 2 a+2 K+2 , i = {a, b}. First, 
we note that for k > and n£N, the following series converge: 

^ (l+|i|)«+i(l + |n-i|)"+3 ' ^(1+ |fc|) K+3 (l + |n-/c|)«+ 1 ' ^ 
Consequently, we have for all (n,k,l) € Z 3 : 

1 IIOn)n6Z ; S a+2 , K+ 2|| ||(9r7n)nSZ ! B a+1 , K+3 \\ 



r 3+2a (l + |||)«+l(l + | n _||)«+3 

and 

\l, a-l \Fl b I \\ s 1 ll(7n)neZ > #a,K+4 1| || (<9 r U>„)neZ ] #a+3,K+l| 
\ K lk\ r )°rW n _ k {r)\ ^ 



r 3 + 2Q (l + |fc|)»«+3(l + | n _jfe|)/e+l 



and therefore the series defining F n is converging. We now bound the series (|33[) . By symmetry, it is 
sufficient to consider only the first series and n > 0. We split the sum into two parts: 



S 



V2] 



n/a (l + m)«+i(l + | n -/|) K + 3 

I — — OO 

<c r 1 v + v 1 <c r 1 x 

- K \l + \n\) Z^(l + |Z|)«+i " K \l + \n\ 
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and 



Q+ = V - 

° n /2 (i + m)«+i(i + | ?1 -m«+3 



<c K 



This shows that F £ 84+20, k+i- 



1 



1 + \n\ 



K + l 



^ a + u 



(1 + H)«+ 3 - Ck v 



5.1 Proof of Lemma [7] 

From now on, wc assume k > 0. Wc recall that: 

r 2 n 1 

Cn := In + i/J,n\ 2 , Vn € Z. 

In what follows, we use without mention the following properties of £ n : 



1+ (n 



Vn eZ\{0}, 



and 



Xm(Cn) 



\n\ 




(- 


V2 


Vn 


n 




(" 


V2 


Vn 



2\ 2 



VneZ\{0}. 



(34) 



We note that is an increasing function of |n| so that its minimal value (over n £ Z \ {0}) is reached 
for n = ±1 and is equal to: 



Let 



:= - min(p M -2,1). 



For /i > [i cr it we have > 2, so that > 0. We choose a £ (0, a^) from now on. This is the smallness 
condition that is mentioned in Lemma [7] 

With the above conventions, we first analyze the equations which determine w: 
Proposition 14 Given F £ 820+4. k and w* £ £>° +2 , the equations: 



w n (r) = —r- + 
w {r) = - 



00 sF n (s) (r\<- f r sF n (s)f S ^ 

OG -1 pQO 



ds . 



with: 



tF (t) dtds, 

sF n (s ) fl^ C 
1 2C„ 



(35) 
(36) 



ds. 



define C 2 functions. Moreover, we have w £ U^ +2 K+2 > an d there exists a constant C a ,^ < 00, depending 
only on a and (i, such that 

II < Co^ [||(F„)„ eZ ; B i+ 2o, K+ 2\\ + HKW ; B° K+2 \\] . (37) 
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Proof. We only prove (|37[) ; existence and continuity follow in a straightforward way. We first 
treat the case n / and then the case n = 0. Throughout the proof, we use the shorthand Mp for 

||(-fn)neZ ! $4+2a,re||- 

Case n ^ 0. We split the expression defining w n according to (|35l) : 

^n« = J+/r(r)+J^(r). 



By definition of the norm on the space £>4+2a. «, we have: 



(1 + H)«S 4 + 2q ' 



Vs > 1 . 



Using that £„ > p M > a + 2, and that c\n\ < |C«| < C(l + v |m|)|«| and |n| < £„ < (1 + v M)M f° r large 
values of |n|, we get that: 



\IUr)\ 



(?) 



< 



2C„ 



ds 



1 



— / s 1+ «"-( 4+2Q )ds 

(l + N)«|Cn| 



(l + r£ n -(2a+2)\ 



(l + H)" 



-2 



(38) 



for all r > 1. Here we used that, by our smallness condition on a, we have £„ — (2a + 2) > — 1. We also 
have: 



2C„ \s) 



ds 



< 



<C a 



(l + H)"|Cn| 

Mp 



a-^-(2a+4) ds 



1 



' M (l + |n|)«+ 2 r 2 «+ 2 ' 
for all r > 1. Using these bounds for r = 1, we obtain: 

|w n | < \w*\ + C a ,n 



(1 + \n\Y+ 2 ' 

Plugging ([3"g]), (|3"9"]l and (gDJ) into (J3S]) and recalling that > a + 2 yields: 



K(r)| < C a 



M F + \\w*;B° K+2 \ 



(l + |n|) re + 2 

Differentiating with respect to r, we obtain: 



K* + 2 



Vr > 1 



d r w n (r) 



^ + C fm-^IM, Vr>l 



(39) 



(40) 



To summarize, when we differentiate w n with respect to r, the decay in r increases by one power, and the 
decay in n decreases by one power. This observation allows us to bound d r w n in the indicated function 
spaces. Finally, since the expression defining w define a solution of (|16p . we plug the bounds on w n and 
d r w n into this equation and get a bound for d rr w n (r). We obtain that there exists a constant C QiAI , 
depending only on a and fx, such that 



r 2 \d rr w n {r)\ r\d r w n (r)\ 

(1 + M) 2 (i + H) 



+ KM I < Ca, A 



Mf 



K ; ^ +2 I 



i 



(i + H) 



K+2 



r a+2 



Vr > 1 
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We emphasize here that the constant C QjAI depends on a and (i. Nevertheless, it is clear from the com- 
putations above that, when a is fixed and fi varies in a a compact interval Jci, this constant remains 
uniformly bounded. 

Case n = 0. Proceeding as in the case n ^ 0, we get the bound: 



oo 



1 f°° 1 , , „ M 



s 



\Mr)\<M F j -J —dtds<C a ^, Vr>l. 



£4+2 



Similarly, one shows 

A/f /*°° 1 A/f 
\drMr)\<-^l ^ds<C a -^, Vr>l, 

and we again conclude, by recalling the differential equation satisfied by wq (see (|16p for n = 0), that: 

2|» CM. Ifl / M , I /M/r- Mj " + |1 i g° +2 ll w . 1 

r \d rr w {r)\ + r\d r w (r)\ + \w {r)\ < C a ^— , Vr > 1 . 

This completes the proof. ■ 

We next consider the equation satisfied by 7: 

Proposition 15 Given <fi G B a +2,n+2 o/nd 7* G ^^+4; ifoe equations: 

7-M=i- r^#f-V'"d»- r^n N d., (4i) 



ri"l 7 r 2|n| \s/ .A 2|n| 

kCXD 1 />00 



7o (r) = / - / tMt)dtds, (42) 



s 

with 



ln = ln+ k "WW ds ' 

define C 2 functions. Moreover, 7 € Z-/^ K+4 and there exists a constant C a < 00, depending only on a, 
such that 

||(7n)n 6 z ; < K+4 || < C Q [||(^)„ e z ; + ||( 7 *)r ie z ; B° +4 ||] . (43) 

The proof is identical to the proof of Proposition [T4l and is left to the reader. Lemma [7] is a straight- 
forward consequence of Proposition [15] and Proposition [14] 



5.2 Proof of Lemma I12L first item 

Let I = C (jUcHt)Oo) and a be given by (piT]) . In particular, we have a < rnin{a! M ,/x € 7} 

so that, applying the results of the preceding section, it follows that 5 M : £? 4+ 2a, K X (£>° +4 x £>£ +2 ) — >• 
k+4 x k+2 i s a well-defined continuous linear map, for all values of (i £ /. We now show 

that the map 5 : /j t— > 5^ is also continuous. This amounts to show that, for arbitrary (jq € I and 
[i G /, there exists a constant C M which converges to zero as fi converges to /Uo, such that, for arbitrary 
(F,i*)e6 4+ 2«x(B» +4 x6» +2 ): 

\\S^ (F,x*) - S^(F,x*) ; „ +4 x U 2 a+W \\ < CJ(F,x*) ; £ 4+2ct , K x (£° +4 x B° +2 )|| . 

Given n G /, we let (7[/i], := S^(F,x*). Since 7 [/it] is obtained from u>[/u] uia an equation which 

does not depend on fi, we can apply directly Proposition 1141 yielding that, for arbitrary (/i, ft) G I 2 : 

IItN - tIm] ; ^a,«+4ll < c a ||«)M - u>[m] ; ^ +2iK+2 || ■ 

Hence, it suffices to prove that w[jj] is continuous with respect to /i, in order to obtain the continuity of 
S. 
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To show the continuity of w\p], we first remark that tool/*] does not depend [i, so that we only detail 
the case n^O. Let n ^ 0, we split io n [/i] into three terms: 



where: 



w n M(r) = W b n [fx](r) + JTM(r) + I^M(r) 
f r sF n {s) /8\CnM 



1 1 2CnM 



sF n {s) fr\CM 
KM 



ds. 



We recall that 



T* r r* 



9 rr u;„[/i](r-) = F„(r) - 



d r w n [n](r) (n 2 + ifin) 



(44) 



(45) 



Note also that CnM = (n 2 + ifj,n) 1 / 2 is a continuous function of uniformly in n. Indeed, since the square 
root is analytic in a neighborhood of 1, we have for sufficiently large n (uniformly in /i £ 7): 



ICnM - Cn[£]| 



1 + V 



III 

1 + — 



We also have the bound ICnMI < c | tt. | , with c independent oi fi £ I. Introducing these bounds into (|4"4"]) - 
(|45[) shows that the continuity of w\p] follows from the continuity of (W&[/z], I\ [p], Ica[lA) in B a +2,K+2- 
For consistency, the three sequences, which are only defined for n^O, are completed by for n = 0. 

To begin with, we consider the continuity of /i i— > Let (/x, /i) <E / 2 , and assume that |Cn[//] — 

Cn[A]| < a /2, uniformly in n. We have: 



7TM-^[A] = Ji(r) + J 2 (r) 



where: 

Jr(r) = 
■Mr) = 

We have, uniformly in n: 

Cn[#] - CnM 



r sF n (s) 

i 2C„M 

gfn(a) 
2Cn[Al 



(?) 



CM 



(;) 



ds , 

r 



ds 



Cn[£] 



<C|C„M-Cn[l"]|=0(l), 



< c|CnM ~ Cn[A]|lnWr |C " M - fn[/I11 = o(l)ln(r)r a / 2 . 

We introduce these uniform bounds in J\ and Ji and redo the computations in the proof of Proposition 
Ql(see (ESI))- We get: 

||fig4+2a|| o(l)ln(r) 
1 U Jl - (1 + H K + 2 ) r 2 + 2Q ' 1 21 )] ~ (1 + |n| K + 2 ) r 2 + 3 "/ 2 ' 

where the term o(l) denotes a constant converging to when /j, — fi — > 0, uniformly in n . Finally, we 
have that, for all n £ Z \ {0}: 



l/ny _ ™, < H^ g 4+2a|| o(l)_ 

hWl ~ (l + \n\-+ 2 ) r 2 +« ' 



(46) 
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We now prove the continuity of /i i— > .Zoo [/■*]• F° r an Y (p-i fi) G ^ an d n / 0, we perform a similar 
splitting: 

CM = Ji(r) + J 2 (r), 



where: 



Ji(r) 
Mr) 



sF n {s) 

r 2 C„M 

sF n (s) 



Cnlfl] ~ CnM 



As in the preceding bound we have, uniformly in n 



, r\("H-(n[(ir 

1- - 



ds . 



ds. 



1 - 



C„M-Cn[A] 



<c|c„M-Cn[A]| = o(i). 

(;)(;)""• 



where we have used that \Cn[n] — Cn[A]l < We can therefore redo the computations in the proof of 
Proposition QH (see (J3SJ). This yields, for all n e Z \ {0}: 



|Ji(r)| + |J 2 (r)|< 



|F;B 4 +2a,«|| o(l) 



(1 + H K + 2 ) r 2 +« ■ 

As in the preceding estimate we conclude that: 

1 1 - -£x>[/i];#a-+2,K+2|| = o(l)||F;S4+2a, K || ■ 

Finally, we prove the continuity of [i i— > Wb \p] : 



< 



< 





" W n [p] 








r 




W n [fj] ' 


- w„ [/I] 




W„[/2] / 






+ 




\Wn\}A " 


" «U/2]| 




w„[/2]| 1 



1 1 



CnlH-] \. r CnM-Cn^-] r C» [A]"C» [**-] 
1 1 



r.2+2o 



p 2+q 



7-CnM-Cn[M-l rfn[P]-Cn|>*-] 



where we have used that 7?.e(Cn[/i]) > 7£e(£„[/t_]) > 2 + 2a. At this point we note that the bound which 
we obtained above for 1^ in r = 1 yields: 



\w n [n] - w n [p]\ = 



0(1) 



(l + \n\Y 



■^\\F;B 4 +2a,K\ 



As /i i— > C«[m] is continuous in /x (uniformly in n) and satishcs 72.e(Cn[A*]) > 7^ e (Cn[M-])j f° r au M G /, wc 
also have: 



1 f 1 1 \ r oo n V 

where o(l) is uniform in n. By combination, this yields, for all n G Z \ {0} 

o(l) 



'(1), 



\wm{r)-wm{r)\< Xl ^ 

This completes the proof of the first item in Lemma [T2 



F-B i+2a ^\\ + \\w*-Bl +2 \ 
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2 - \n\ - („ 



5.3 Proof of Lemma I12L second item 

In this paragraph, we prove that the map 71 is one-to-one and onto with a continuous inverse. Given 
x* = (ff*,w*) € £>" +4 x B® +2 , we set (7,«)) = 5^(0, x*). A straightforward computation shows: 

f , K 7» , f°° sw n (s) /r\M /"" sw„(s) /s\ H 

where: 

Therefore, we have, for all n e Z \ {0}: 

7n(l) = 7^ » 

together with: 

f°° /1\ H 
9r7n(l) = ~Wnl + J sw n {s)\-\ ds = -\n\j* 

so that (v*,Vg) =Ti(x*) satisfies: 

• K,o = v *e,o = 

• v ; n = inj*, and v* 6 n = -d rln {l) = |n| 7 * + - _ J_ - , for all n € Z \ {0}. 

This shows that the map 71 is one-to-one and onto. Indeed, the inverse map is given by: 

T 1 - X [v' r ,v' e )] = {i,w'), 

where: 

• To = w 'o = °> 

. i n = and w' n = (2 - |n| - Cn) U n - for aU » G Z\ {0}. 

It is therefore clear that T^ 1 € £ c (£J£ +3 x B° +3 ; S° +4 x B° K+2 ). This completes the proof. 
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